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'sj" . Abstract 

in 

I We present a simple interpretation of the recently observed near- 

■ threshold proton- antiproton enhancement. It is described by a set of 

T^lj- , low-energy parameters deduced from the analysis of NN experiments 

I at LEAR. We predict a related effect in photoproduction reaction 

• under study by CLAS collaboration. 

{iLiI Low-mass baryon-antibaryon enhancement has recently been observed in 

^ ; the decays B+ K+pp [I], 5° ^ D^pp [2] and J/ip 7pp 0- On the 

theoretical side this discovery has been discussed by several authors [U El El 
. E|. In 13] the near-threshold effect has been considered either as a gluonic 

! state or as a result of the quark fragmentation process. In [Hj it has been 

I attributed to the Breit colorspin interaction while in jB] it has been regarded 

as caused by peripheral one-pion-exchange potential. The quantum numbers 
corresponding to the observed enhancement are discussed in 

The purpose of the present work is to show that the recently observed 
near-threshold pp enchancement may be understood invoking the knowledge 
on NN interaction gained at the Low Energy Antiproton Ring (LEAR) 
at CERN. The LEAR results may be schematically subdivided into three 
groups: 

(i) meson spectroscopy, 

(ii) baryonium searches, 

(iii) NN interaction: scattering, annihilation, protonium spectroscopy. 
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A comprehensive review of the investigations performed at LEAR along 
these hnes may be found e.g. in [S]. At first sight it might seem tempting 
to connect the newly observed structure with the elusive baryonium [9 -^JJ. 
The point of view adopted in this paper is different. We remind that bary- 
onium searches at LEAR ended without any clear evidence for its existence 
[Hj. On the other hand the studies of NN elastic scattering, annihilation and 
charge-exchange reaction along with the protonium studies yielded reliable 
information on low-energy NN amplitude [5*. It has been shown [T^-[T3] 
that the whole set of the experimental data on low-energy NN interaction 
is possible to describe in terms of the effective range approximation. In 
this approach the lack of the dynamical picture is traded for the possibility 
to reconcile within a unique scheme different pieces of information on NN 
interaction. 

We shall see that the effective range analysis based on the old 

LEAR data enables to explain Belle ^1 12] and BES |3] results. In particular 
the observation by BES of a strong effect in the decay J/ip ^ jpp and the 
absence of a similar structure in J/tp ^ n^pp perfectly fits into the solution 
for the low-energy parameters obtained 15 years ago ^21-^2] • We shall return 
to this point below. 

The problem of whether any conclusions on baryonium can be inferred 
from the low-energy parameters is beyond the scope of this work. We only 
mention that the values of the scattering lengths presented below are not 
anomalously large as it should be for the case of a level close to threshold. 
One must also keep in mind that strong annihilation drastically changes the 
standard relation between the low-energy parameters and the positions of 
the poles of the amplitude p!B] . 

Our approach is based on the low-energy analysis of final-state pp in- 
teraction. Therefore the near-threshold enhancement observed in different 
reactions jU |2| is described by the same equations presented below. Another 
consequence of the proposed scheme is the prediction of the similar phenom- 
ena in other reactions of the same type, in particular in photoproduction 

— > ppp. This process can be investigated by CLAS collaboration at Jef- 
ferson Lab. Unlike B- and J/ip decays the photoproduction reaction has not 
been discussed from the theoretical side. Therefore we shall choose it as a 
starting point in order to introduce the effective range formalism. Then we 
shall return to Belle and BES data. 

The double differential cross section for the reaction jp ppp is given 
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by the well-known Chew-Low expression 



ds2dt 
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where m is the nucleon mass, k is the energy of 7, 



(2) 



with P2 and ps being the 4-momenta of p and p forming the low-mass pair, pi 
being the 4-momentum of the remaining proton, \{x, y, z) = {x—y — zY—'iyz 
is the standard kinematical function. The angle (^23 = {cos9,ip) is defined 
in the CM system of the particles (2,3). The variable S2 = Tn\^ is the square 
of the invariant mass of the (2,3) pair. It follows from (jS)) that the angular 
distribution of particles 2 and 3 ( the low-mass pp pair) is isotropic in their 
CM system provided T(s2, ti, ^^23) independent of (^23 which is the case 
for ^5*0 and ^Pq states. Equation (jT)) yields the following invariant mass 
distribution 



where (3 is the angle between and the direction of the incident 7 ( we 
remind that the index 1 is attributed to the proton which remains outside 
the correlated pp pair). 

In the spirit of the Migdal- Watson FSI theory we single out from the 
matrix element T a factor responsible for the low energy pp interaction. 
Within Migdal- Watson approach one has 



where q is the pp CM momenta, = ^{ml^—Am'^), f{—q) is the Jost function 
corresponding to the pp interaction at low energy, D{q) = |/(— g)]^^ is called 
the enhancement factor. As it was already mentioned the dynamics of NN 
interaction is much more complicated than that of NN jS], [0], ^I]; [TH] . 
Annihilation dominates at short and possibly intermediate distances, u- and 
other odd G-parity exchanges lead to a strong attraction to be added to 
the one-pion exchange in the outer region 117]. Therefore the approximation 
of f{—q) by a Born term from one-pion exchange jE] misses the essential 
features of NN dynamics. On the other hand effective range solution for 




(3) 



T|2 = |T(0)|V|/(-g)P = D(g)|r(°)p, 



(4) 
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the low-energy NN amplitude describes the whole set of data with a fair 
accuracy [121- [13 • 

Ignoring for the moment complications due to annihilation, Coulomb in- 
teraction and spin-isospin structure we may write the following expression 
for Jost function in the scattering length approximation 

f{-q)^A{q){l-tqa), (5) 

where a is the scattering length and A{q) is a well defined smooth function. 
A similar expression for pp system is not as simple for three reasons: 

(i) pp system is a combination of the two isospin states with I = 0, 1; 
isospin invariance is violated by the mass difference of proton and neutron; 

(ii) a powerful annihilation results in the complexity of the scattering 
lengths; 

(iii) the Coulomb attraction acts in the pp system. 

Experimental data on NN spin observables exist only at the incident 
p momenta Pl > 500 MeV/c fH^ IHj. i.e. beyond the energy region where 
the effective range approximation is applicable. From the level shifts in pp 
atom it is possible to extract only some information on the imaginary 
parts of the triplet and singlet scattering lengths ^Uj. Therefore only spin- 
averaged effective range parameters were extracted from the experimental 
data [I2]-[in]- 

The enhancement factor in which all three points (i) -(iii) inherent for the 
pp system are accounted for has the following form |12j-[T3] 

^^^^ " \l - ts{q + lA + I) - (g + zA)/r|2' 



where 

c (g) 





27T Y 


1 — exp ^ 






qaBj _ 



. . (7) 

gas 

is the Sakharov Coulomb enhancement factor |2I] with = 2/ am = 57.6 
fm being the Bohr radius of the pp atom, 

■5 = ^(«o + «i), r = aoai, (8) 

where aj, 7 = 0,1 are NN S -wave scattering lengths with isospin /, 

q = c^(q)q+ — hiqan), h{z) = \n z + Retp (l + - ) , (9) 
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with ip{z) = d/ dz In r(2;). The quantity / is the momentum in the nn channel, 

f = q^-m5, (5 = 2(m„-m), (10) 

(recall that m is the proton mass). 

The point q = (m5)^/^ ~ 49 MeV/c corresponds to the nn threshold; 
below this point nn momentum / becomes imaginary, / = i(mS — q^Y^'^. 
Parameter A is the Schwinger correction to the scattering length [22], A ~ 
-0.08 fm [I2]-[IS1. 

The enhancement factor ® does not include the effective range term 
and the contribution from nonzero orbital momenta. According to [T^-[T3] 
the effective range term is of minor importance up to g ~ 150 MeV/c i.e. 
io Q = — 2m ^ 20 MeV. We note in passing that for a multichannel 
system the effective range may be even negative or complex (221 121|- The 
P-wave comes into play a little earlier especially as far as the total cross 
section is concerned. Inclusion of P-wave brings about two problems. First, 
the P-wave scattering length (it has a dimension of fm^) is very sensitive to 
the fitting procedure (see |I2]-|ISI) and the present set of the experimental 
data do not warrant a stable solution for the / = 1 amplitude. The second 
problem is the following. With P-wave included equations become rather 
cumbersome and the number of parameters increases substantially. This 
may cause unnecessary doubts in the reliability of the proposed approach. 
We plan to consider effective range terms and / > amplitudes in the next 
publication. This would allow to analyze the pp correlation function in a 
wider energy range. 

The most remarkable result of all fits 52]- CH] for low energy NN param- 
eters is a clear dominance of the S'-wave with / = over that with 1 = 1. 
The corresponding scattering lengths read 

ao = (-1.2 + i0.9) fm, ai = (-0.1 + iO.4) fm, (11) 

ao = (-1.1 + Z0.4) fm, ai = (0.3 + iO.8) fm, (12) 

where (fTT|) and (fT^ are respectively the results of [13] and [U]. The sign 
convention in ll2j-||TH] is A; cot 5 = +l/a. We see that the absolute value of 
|Peao| is much larger than that of |Peai|. This is completely in line with 
the observation by BES of a strong low-mass effect in the decay process 
J /ip ^ '-fpp and the absence of a similar structure in J/ip Tc^pp [3J. 

The two solutions given by (fTT|) and (fT^ are somewhat different. This is 
due to several reasons described in p^-||15j. The main source of ambiguity 
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is the lack of the experimental data at very low energies. The effective 
range analysis which resulted in the parameter sets (|ll |) - (jl2|) was performed 
at the end of eighties. In view of the LEAR shut down the experimental 
situation did not drastically change since that time. The only substantial 
new piece of information concerns the total and annihilation pp and np cross 
sections down to pLab = 35 MeV/c [2nil2Sl- However in the analysis of these 
data [23 I2H] use was made of the effective range solution [22] performed in 
1988 which was later updated in 1991 ^1] partly by the original authors. 
The parameters set (fT^ is just this updated version. The effective range 
solution fairly well fits the low energy cross sections [23]-j2H] which is 
an indication that the same is true for our parameter set (|12p. We urge our 
colleagues who are in possession of the experimental data to persuade the 
analysis whith the parameter sets fT^.T^. An independent phenomenological 
analysis of the experimental results ^25, was performed in [30j. The authors 
extracted the imaginary part of the Coulomb-distorted pp scattering length. 
This imaginary part corresponds to the averaging over spin and isospin. The 
nn threshold which lies just within the region under investigation was not 
taken into account. This fact should be kept in mind when comparing our 
parameter sets which the value presented in jSO]. The result of is Imacs = 
0.69 fm while from ()11|) we get |/m(ao+ai) = 0.65 fm (the index in a^s means 
that this quantity still contains the Coulomb correction) . Again we conclude 
that the "new" data [23 do not contradict the "old" solution (fTTj). 

The negative sign of Rea^ may be interpreted as an indication that the 
potential in this channel is either repulsive or on the contrary attractive and 
strong enough to produce a bound state somewhere below the threshold. 
However such a simple interpretation is not obvious in strong annihilation 
regime. This problem is beyond the scope of the present work. 

The above arguments based on the analysis of the reaction 7p — > ppp led 
us to the following conclusion. Within the framework of the FSI theory and 
within the region of the applicability of the scattering length approximation 
the pp effective mass distribution is governed by the function 



nQ) = (ml 



pp 



AmY^^DiQ) ^ 2{mQ) 



(13) 



where 



Q = m. 



'PP 



2m 



2m 




.2 \ 1/2 



1 



(14) 
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and D{Q) is defined by The phase-space factor Q^/^ is trivial and plays 
the role of the universal background. The function D{Q) reflects the the 
main features of the pp dynamics. Due to Coulomb attraction it diverges at 
threshold since C^(g) — >■ 27i/aBq as q —>■ 0. Thus at Q = the distribution 
2{mQY^'^D{Q) starts from a finite value An /as- At q = [mSy^'^ the factor 
D{Q) has a cusp due to the opening of the nn channel. 

The above conclusions also hold true for the reactions —>■ K~^pp^, 

_^ d^ppI^ and J/ip 7j9p[3j. 

An remark is due at this point. The BES data on J/ip 'jpp were 
analyzed making use of the Breit- Winger formula [3] . We wish to remind that 
in the vicinity of the threshold the Breit- Wigner function can be reexpressed 
in terms of the effective range parameters |31j . 

Our main results are presented in Figs. 1-2. In Fig. 1 we display the 
enhancement factor D{Q) corresponding to the parameters set (fTT|) . In this 
figure we also plot the experimental points obtained by BES collaboration 
and taken from Ref.p]. The normalization of the experimental points is 
chosen in such a way that in the region of Q ~ 0.02 GeV they are in rough 
agreement with the theoretical curves. The comparison therefore may serve 
only for the orientation purposes and in order to give an impetus for the 
future detailed analysis of Belle and CLAS data as soon as these data are 
available. 

The curve D{Q) corresponding to the parameter set (12) is not shown 
since it goes very close to the presented in Fig.l. In Fig. 2 we illustrate the 
roles played by the Coulomb attraction and annihilation. The solid curve in 
Fig. 2 is the same as in Fig.l, the dashed one corresponds to the Coulomb 
interaction switched off, while the dotted curve - to switching off the anni- 
hilation. 

Our main conclusions may be summarized as follows. The near-threshold 
proton-antiproton enhancement has been observed in reactions driven by 
different mechanisms. Therefore the effect may be due to pp FSI. In order to 
verify this assumption we have described the FSI in terms of the low-energy 
NN parameters taking into account the isospin structure of the amplitude 
and the Coulomb interaction. The factorization of the FSI factor from the 
total amplitude has been explicitly performed for the reaction jp — >■ ppp. In 
turns out that the data of BES collaboration [3] are fairly well described by 
the effective range solution obtained long ago from the analysis of the LEAR 
data. Further conclusions would become possible as soon as new experiment 
data from BES, Belle and CLAS are available. Finally we remark that similar 
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kind of analysis may be applied to AA final states. The low energy parameters 
for this system may be found e.g. in [22]. 

The authors are indebted to L.N.Bogdanova, A.Kudryavtsev, K.Mikhailov, 
P.Pakhlov, Yu.Simonov and B.S.Zou for useful discussions and clarifying re- 
marks. We are grateful to Dr. Shi-Lin for drawing our attantion to paper 
(7j. B.K. acknowledge the support from the grant Ssc-1774-2003. 
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Fig. 1 The near-threshold behaviour of the enhancement function D{Q), Q = m23 — 2m. 
Sohd curve corresponds to the parameter set (fTTj) . Experimental points are from Ref. [3]. 
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Fig. 2 The enhancement factor according to (jHl) with parameter set - the sohd curve, the 
same without Coulomb interaction - the dashed curve, without annihilation - the dotted curve. 
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